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We study timelike and spacelike virtual Compton scattering in the generalized Bjorken scaling 
regime at next to leading order in the strong coupling constant, in the medium energy range which 
will be studied intensely at JLabl2 and in the COMPASS-II experiment at CERN. We show that the 
Born amplitudes get sizeable 0(a s ) corrections and, even at moderate energies, the gluonic contribu- 
tions are by no means negligible. We stress that the timelike and spacelike cases are complementary 
and that their difference deserves much special attention. 

PACS numbers: 12.38.Bx, 12.38.-t, 13.60.Fz 



INTRODUCTION 



C~| ' Spacelike Deeply Virtual Compton Scattering (DVCS) [l[ 

1 *(q ln )N(P)^ 1 (q out )N'(P / = P + A), qf n = -Q 2 , q 2 out = , (1) 

has been the model reaction for studying QCD collinear factorization in exclusive processes in terms of Generalized 
i . Parton Distributions (GPDs) 0, H| , which access correlated information about the light cone momentum fraction and 
transverse location [4| of partons in hadrons. The specific kinematical regime where this factorization property is 
proven at leading twist is the generalized Bjorken regime of large energy and large Q 2 but finite and fixed momentum 
transfer A 2 . A number of experimental results at various energies [B| have now established the relevance of this 
approach at accessible kinematical conditions. Detailed phenomenological studies @ are under way to quantify to 
which degree one may in a foreseeable future extract from experimental data the physical information encoded in 
GPDs. The DVCS process contributes to the leptoproduction of a real photon 

^(fe)iV(P) -y ^(fc'h^onO W = P + A) ■ (2) 
d | It interferes with the Bethe-Heitler process 

> : l ± (fc) 7 *(-A) -> ^(fc'h^t) , (3) 

• rH ' 

where the hadronic interaction is entirely determined by the nucleon (spacelike) electromagnetic form factors F\{A 2 ) 
^ andf 2 (A 2 ). 

Timelike Compton Scattering (TCS) Q 

1 (q m )N(P)^ 1 *(q out )N'(P' = P + A), q 2 n =0, q 2 out = Q 2 , (4) 

which contributes to the photoproduction of a lepton pair 

1 (q m )N(P)^l-(k)l+(k')N'(P' = P + A), k + k' = q out , (5) 

and interferes with the Bethe-Heitler process 

7fe„) 7 *(-A)^r(fc)/ + (fc') , (6) 

shares many features with DVCS and allows in principle to access the same GPDs. The experimental situation [8] is 
not as encouraging as in the DVCS case but progress is expected in the next few years. The amplitudes of these two 
reactions are related at Born order by a simple complex conjugation but they significantly differ at next to leading 
order (NLO) in the strong coupling constant a s Q. 

Complete NLO calculations |10l - fl4l ] are available for both the DVCS and TCS reactions and there is no indication 
that NLO corrections are negligible in the kinematics relevant for current or near future experiments. 
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In this paper, we explore the consequence of including NLO gluon coefficient functions and NLO corrections to the 
quark coefficient functions entering the DVCS and TCS amplitudes, firstly in the calculation of spacelike and timelike 
Compton form factors with two models of GPDs. Then we proceed to the calculation of specific observables in the 
kincmatical conditions which will soon be accessible in lepton nucleon collisions. 

We focus on kinematics relevant to the next JLab and COMPASS measurements. The case of very large energy, 
and subsequently very small skewness, deserves its own study that will be addressed elsewhere. Ultra peri phcral 
collisions at hadron colliders may already open the access to timelike Compton scattering in this domain [15| . Also 
the proposed next generation of electron-ion colliders [161 Il7j j will access in a detailed way this interesting domain 
with small skewness. 

In this paper we concentrate on the influence of the NLO corrections to the DVCS and TCS observables. Because 
of that, although required to ensure QED gauge invariance [l8j], we do not discuss here twist 3 effects. Neither do 
we implement the new results on target mass and finite t corrections [l9j nor the recently proposed resummation 
formula [2(| which focuses on the regions near x = ±£. 



2. KINEMATICS AND AMPLITUDES 



1. Kinematics 



We introduce two light-like vectors p and n satisfying p 2 = 0, n 2 = and np = 1. We decompose the momenta 
in this light-cone basis as k^ = an^ + (3p^ + fcj., with k\~, < 0, and we further note k + = kn. We also introduce the 
standard kinematical variables: A = P' — P, t = A 2 and W 2 = ((&„ + P) 2 . 

We define kinematics separately for spacelike (si) and timelike (tl) Compton scattering. We denote the (positive) 
skewness variable as £ in the DVCS case, and as r\ in the TCS case. In the DVCS case, where — qf n = Q 2 > and 
lout = 0, we parametrize the momenta as follows 

<4 = ^ - ^P" , lout = a s m» - - A£ , 

M 2 M 2 — A 2 

P» = (l + ey + — — ^, P**=/W+ 0R T n^ + A£, (7) 
2(1 + £) Zpsi 

where M is the nucleon mass. The coefficients a s i, f3 s i in satisfy the following system of equations 

Q 2 M 2 ( 1 1 \ A 2 A 2 

a sl = -jr H r- ( — -t - -r + jrf- , Psl = 1 - £ + ^ , (8) 



which in the limit M — and Ay = relevant for calculation of the coefficient function leads to the standard values 
a s l = Q 2 /4£ and p ai = l-£. 

In the TCS case, where q 2 n = and q 2 out = Q 2 > 0, we parametrize the momenta as 



O 2 O 2 - A 2 



M 2 M 2 - A 2 

P" = (1 + r/K + -— n" , P'^ = ft^ + - T n» + A£ . (9) 
2(1 + ?7) 2/3 t ; 

The coefficients au , /?tz in ([9]) are solutions of the following system of equations 

Q 2 M 2 ( 1 1 \ A 2 Q 2 - A 2 



4ry 2 Vl + ^ PuJ Wti ' " ' 2a t 

which, again in the limit M — and Ay = relevant for calculation of the coefficient function, take the standard 
values au = Q 2 /^"q and flu = 1 — V- 
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2. The DVCS and TCS amplitudes 



After proper renormalization, the full Compton scattering amplitude 1 reads in its factorized form (at factorization 
scale up) 



-gt^" J dx 



ICr, 



dx 



J2T q {x)F q (x) + T 9 (x)F 9 (x 

Q 

l_F ^ 

^T q {x)F q {x)+T 9 {x)F 9 (x) 



(11) 



where we omitted the explicit skewness dependence. Renormalized coefficient functions are given by 

12' 



T q {x) 
T q (x) 
f q (x) 
f q (x) 



C q {x) + C q {x)+\n(^-C q coll {x) 



{x — > —x) , 



C 9 {x) + hx i9L\.ci o Jx) 



+ (x — >■ —x) 



C q (x) + C q (x)+ln )-C q oll (x) 



Q 2 



+ (x —x) , 



C 9 (x)+\n[^-)-C 9 oll (x) 



— (x — > —x) 



(12) 



Results of the NLO calculations 

[IM1 of the quark coefficient functions, based on the standard definitions of GPDs 
given in the Diehl's review 0, read in the DVCS case 



C?(x,0 
C q (x,0 





1 

IE 




e 2 q a s C F 


1 




47T X - 




IE 


e 2 q asC F 


1 




47T X - 


-e- 


IE 


e 2 1 






€q x + t- 


IE 




e\a s C F 


1 




47T X - 




IE 


e\a s C F 


1 




47T X - 




IE 



9 - 3 log(^-— - is) - log ( 

x-£ 2£ 



3_21og(^-z £ ) 



2£ 



9 -^ log( ^r 

_3-21og(H±l-fe) 



(13) 



where CV = (TV 2 — l)/(2iV c ), N c being the number of colors, and e q is the quark electric charge in units of the proton 
charge. Using the same conventions, gluon coefficient functions read in the DVCS case 



Cf(x,0 = 



T,e 2 q asTF 



1 



47r (x + £ — ie)(x — £ + is) 



X + 3£. fx + t; 



T,e 2 q a s T F 



log 



2£ 



IE 



X + Z , 2 ( X + Z 



x-i 



lo. 



2e 



x + /x + £ 



47r (a; + £ — ze)(x — £ + ie) [ x — £ 
£e 2 asT F 1 

47r (x + £ — ie)(x — £ + is) 



lo. 



2£ 



1 We do not consider the photon helicity changing amplitude coming from the transversity gluon GPD . 
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- 2 log — — l£ H log — — IE 



=2 

9 ^ ^ - rr« 



(14) 



47r (a; + £ — ie)(.T — £ + ie) 
where T F = |. The results for the TCS case are simply Q related to these expressions 

TCS T(x,r,) = ±( DVCS T(x ) Z = r,)+inC cM (x^ = r ] ))* , (15) 
where + (— ) sign corresponds to vector (axial) case. 

3. MODELS FOR GPDS 

In our analysis we will use two GPD models based on Double Distributions (DDs) 0, [22[ • DDs allow to trivially 
achieve one of the strongest constraints on GPDs : the polynomiality of the Mellin moments of GPDs. They also 
automatically restore usual PDFs in the forward limit at £, t — > 0. The GPDs are expressed as a two-dimensional 
integral over a and j3 of the double distribution fa 

Fi(x,Z,t)= [ d/3 f m da6(p + ta-z)fi(P,a,t)+D? (*,t) 9(£ 2 - a; 2 ) , (16) 
J-i J-i+\p\ 



where F = H, E, H , E and i denotes the flavor (val for valence quarks, sea for sea quarks and g for gluons). In our 
analysis we only take into account the contribution of H and H . Indeed, E and E are mostly unknown, and recent 
phenomenological studies of Ref. [23| showed that most of existing DVCS observables are sensitive mostly to H and 
H. 

The DD fi reads 

t(R (RAh(ff\ T ^n % + 2) [(l-|/3|) 2 -a 2 ]"- 

where T is the gamma function, rii is set to 1 for valence quarks, 2 for sea quarks and gluons. gi(f3,t) parametrizes 
the i-dependence of GPDs, the /ij(/?)'s in case of GPDs H and H denote their forward limit and are related to the 
usual polarized and unpolarized PDFs in the following way 

h g W) = \P\g(\P\), 
h q sc M = <Zsca(|/3|)sign(/3), 

K*(P) = 9vai03)e(/?), 

h g ((3) = pAg(\p\), 
hlM = Aq sca (\l3\), 

h%((3) = Ag val (/3)9(/3). (18) 

Df in Eq. (|16l) denotes the Polyakov- Weiss D-term [24j . In our estimates we will use parametrizations obtained by a 
fit to the chiral soliton model [25j | : 

D»q,t) = -D*q,t) = ±D*q)F D (t), 

where 

oo(odd) 

D«(x,t) = -(l-X 2 ) < C h*), 

71=1 

oo(odd) 

D«{x,t) = --(1-.t 2 ) 2 J2 dB n C Li(x)> (2°) 

Z n=l 

where at fi = 0.6 GeV matching to chiral soliton model gives: d\ = 4.0, = 1.2, d\ = 0.4, and we make an assumption 
that at input scale vanishes. In the QCD evolution of d^d\ we switch from 3 to 4 flavours at /i = 1.5 GeV and 
A 3 = 0.232 GeV, A 4 = 0.200 GeV. 
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1. The Goloskokov-Kroll model for the GPDs 

As described in details in Refs. |23l l26H28j | ■ the GPDs of the so-called Goloskokov-Kroll (or GK) model is constructed 
using CTEQ6m PDFs 129| . The low- a; behavior of PDFs is well reproduced by power-law, with the power assumed 
to be generated by Regge poles. In this GPD model, the Regge behavior with linear trajectories of the DD encoded 
in the function gi((3,t) (see Eq. (ITT)) ') is assumed 

g i {S3,t)=e h * t \l3\-< t . (21) 

For the unpolarized GPD H(x,£,t) the values of the Regge trajectory slopes and residues a[,bi as well as the 
expansions of the CTEQ6m PDFs [2!| used for hi may be found in Refs. [26U28j | . Finally, this model uses simple 
relations to parameterize the quark sea 

ttU rrd ttS 

^sea ^sea "'s^sea ' 

with k s = 1 + 0.68/(1 + 0.52 \nQ 2 /Ql), (22) 

with the initial scale of the CTEQ6m PDFs Ql = 4 GeV 2 . 

Similarly to the GPD H, polarized GPD H is constructed using the Bliimlcin - Bottcher (BB) polarized PDF 
parametrization [33] to fix the forward limit. Meson electroproduction data from HERA and HERMES have been 
considered to fix parameters for this GPD in the GK model, H for valence quarks and gluons have been parametrized, 
however we have neglected H for sea quarks. The values of the Regge trajectory slopes and residues as well as the 
expansions of the BB PDFs used in the GK model may be found in Refs. [26|, [31j ■ 



2. The MSTW08 based model with factorized t - dependence 

For the second model we use double distribution with MSTW08 PDFs 32]. In that case we take simple factorizing 
ansatz for t - dependence 

S«GM) = (23) 

g d (j3,t) = F?(t), (24) 

g 8 (fi,t) = g g (P,t) = F D (t) , (25) 

where 

f?(t) = 2Ff{t)+F?(t), (26) 

F*(t) = Ff(t) + 2F 1 n (t), (27) 

F D (t) = (l-t/M*)- 2 , (28) 

with My = 0.84 GeV, Ff and F™ are electromagnetic Dirac spacelike form factors of the proton and neutron. We 
use that model to construct only %. 

4. COMPTON FORM FACTORS 

Let us now present the results for spacelike and timelike Compton Form Factors (CFF) at NLO, T-L and H, defined 
in the DVCS case as 



H&t) 



J \ dx {ll T<? ^ ^ H "^ ^ f ) + T9 ( x ' ^ H9 ^ ^ *)) 
J dx \ Y t T'{x, 0H«(x, C, t) + V>(x, Z)H g (x, £, t) 



H(£, t) = -Jte\^ OH^x, £, t) + TO(x, 0H°(x, £, t) j . (29) 

These CFFs are the GPD dependent quantities which enter the amplitudes. For DVCS they are defined through 
relations such as Q 

A^, t) = -e 2 {p ^ p/)+ u(P') [ (U(t t) 7 + + t) % ^ £ ) + < t) 7 + 7s + t) ^f) ] u(P) , 

(30) 
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Similar relation holds for TCS with £ replaced by rj. 

We now present our results for the spacelike and timelike virtual Compton scattering, for the £ and 77 values which 
include kinematical regimes of the JLab and Compass experiments. As we present our results for Q 2 = = n 2 R = 4 
GeV 2 , throughout the whole paper we use the value of as — 0.3 . 

1. Spacelike Compton Form Factors 

Let us first discuss the importance of including NLO effects in CFFs related to DVCS observables. We show on 
Fig. Q] (resp. Fig. [2]) the results of our calculations for the real (resp. imaginary) parts of the CFF (multiplied 
by £ for a better legibility of the figure, since the CFFs are roughly proportional to l/£), with the two GPD models 
described in Sect. [3] The dotted curves are the LO results and the solid lines show the results including all NLO 
effects. Although the results are naturally dependent of the choice of model GPDs, the main conclusions are quite 
universal. One can first observe that NLO corrections are by no means small, as exemplified by the ratio of these 
NLO corrections to the LO result shown in the lower part of Figs. Q] and [H One can also observe that the NLO 
corrections tend to diminish the real part of the CFF, and even change its sign for £ > 0.01). NLO corrections also 
decrease the imaginary part of the CFF. These are not new results [TlT - [l3| |. To quantify the main source of the NLO 
contribution, we show on the same plots with dashed lines the real and the imaginary parts of the CFF including 
quark NLO effects but not the gluon effects. Gluon effects are the most important part of the NLO correction, even 
at quite large values of £ (up to around £ s=s 0.3) and they contribute as a very significant part to the full CFF "H(£) 
including NLO effects. Since the CFF H(£,) dominates the DVCS amplitude, this means that extracting quark GPDs 
from a leading order analysis of DVCS data is, to say the least, questionable. More positively, this result indicates 
that DVCS experiments even in the low energy regime of JLabl2 provide us with a nice way to measure gluon GPDs. 
This fact has, to our knowledge, never been clearly spelled out before. 

For completeness, we show on Fig.[3]the real and imaginary parts of the CFF £%(£) only for KG model. As for the 
case of %(£), although the effects are less dramatic here, the NLO corrections are by no means small. 

Let us now comment on the D-term contribution to the CFFs. Since the GPD originating from a D-term is a 
function of the ratio x/^ with a support in the ERBL region — £ < x < £, it results in a constant real contribution 
to CFFs in the spacelike case. For Q 2 = /ip = 4 GeV 2 and t — —0.1 GeV 2 , the values of different contributions to 
Ht> are shown in the second column of the Tab. HI One can see that this D-term contribution is very important for 
£ ~ 0.1, since it modifies the the full NLO result by about 50%. It is much less important at lower values of £ (around 
3% for £ - 0.01 and only a few per mil at £ - 0.001). 

2. Timelike Compton Form Factors 

Let us now discuss the corresponding results for the timelike CFFs. On Fig. [4] and Fig. O we show the real and 
imaginary parts of the CFF H for the KG and the MSTW08 models of GPDs described in Sect. [31 for the invariant 
mass of the lepton pair Q' 2 = 4 GeV 2 , t = — .lGcV 2 and factorization scale fip = Q. For the imaginary part the 
correction does not exceed 40%. In the real part, the correction is of the order of few hundred percent. We observe 
that the main part of that large correction comes from the contribution of gluonic GPDs. To quantify the sensitivity 
of this statement on the uncertainties of the input PDF parametrizations we replot on the Fig. [7] the upper right 
panel of Fig. 2] and the upper right panel of Fig. [SJ now with the shaded bands showing the effect of a one sigma 
uncertainty of the input MSTW08 fit on the full NLO result, compared to the LO. 

The D-term contribution to the CFF is a /^-independent quantity and it has both a real and an imaginary parts 
at NLO in the TCS case. We show in Tab. U the values of this D-term contribution in the LO and NLO cases. Its 
relative effect on the imaginary part of the CFF decreases significantly when r\ decreases, from 10 to 1 and 0.1% when 
77 decreases from 0.1 to 0.01 and to 0.001. 

We then compare TCS and DVCS by plotting the ratio of NLO corrections on Fig. [5J There is a striking difference 
in the magnitude of the corrections to the real part of CFFs, mostly insensitive to the choice of GPD parametrizations. 
As discussed in Ref. Q , this is a consequence of Eq. (fT5|) which by adding a phase to the dominant imaginary part 
of the spacelike CFF at small skewness, gives rise to a sizeable real part of the corresponding CFF in the timelike 
case. Such large corrections to the real part of CFFs will have a significant influence on observables which depend 
on the interference of the TCS process with the Bethe-Heitler amplitude, i.e. connected to the azimuthal angular 
distribution of the leptons. We shall discuss this in the next section. 

For completeness, we show on Fig. [8] the real and imaginary parts of 77% (77). The NLO corrections are here smaller 
than 20% in the whole domain. 
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1.0 n ' i ' ' ' ' ■ ■ i ■ ■ " i 1.0 rr 




Figure 1: The real part of the spacelike Compton Form Factor H(£) multiplied by f, as a function of £ in the double distribution 
model based on Kroll-Goloskokov (upper left) and MSTW08 (upper right) parametrizations, for up — Q 2 = 4 GeV 2 and 
t = —0.1 GeV 2 , at the Born order (dotted line), including the NLO quark corrections (dashed line) and including both quark 
and gluon NLO corrections (solid line). Below the ratios of the NLO correction to LO result in the corresponding models. 



2.0 n ■ — ■ ■ j 2.0 n ' — — i— 

1.5 - 1.5 




Figure 2: The imaginary part of the spacelike Compton Form Factor H(£) multiplied by £, as a function of f in the double 
distribution model based on Kroll-Goloskokov (upper left) and MSTW08 (upper right) parametrizations, for n 2 F = Q 2 — 4 GeV 2 
and t = —0.1 GeV 2 , at the Born order (dotted line), including the NLO quark corrections (dashed line) and including both 
quark and gluon NLO corrections (solid line). Below the ratios of the NLO correction to LO result in the corresponding models. 
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Figure 3: The real (left) and imaginary(right) parts of the spacelike Compton Form Factor W multiplied by £ as a function of £ 
in the double distribution model based on KG parametrizations, for = Q 2 = 4 GeV 2 and t = —0.1 GeV 2 , at the Born order 
(dotted line), including the NLO quark corrections (dashed line) and including both quark and gluon NLO corrections (solid 
line). The plots in the lower part show the corresponding ratios of the NLO correction to the LO results. 





ReHv 




LO 


-2.59 





NLO quark contribution 


-0.16 


-0.85 


NLO gluon contribution 


0.18 


0.16 


Full NLO 


-2.57 


-0.69 



Table I: Different contributions to the D-term. The values of the real part coincides for spacelike and timelike CFF H, while 
the imaginary part is non-vanishing only for the timelike case. 



5. CROSS SECTIONS AND ASYMMETRIES 



1. Deeply virtual Compton scattering 

Let us first briefly review the effects of including NLO corrections on the DVCS observables. On Fig. [9] we show the 
total DVCS cross section, the difference of cross sections for opposite lcpton helicities and the beam spin asymmetry 
A^jj, defined by Eq. (48) of Ref. |23|. We choose a set of values of kinematic variables representative for JLabl2, 
namely Q 2 = 4 GeV 2 , E e = 11 GeV and t = —0.2 GeV 2 . The observables are shown as a function of the azimuthal 
angle <j> (in the Trento convention). The Born order result is shown as the dotted line, the full NLO result by the 
solid line and the NLO result without the gluonic contribution as the dashed line. We see that the effects of the 
NLO corrections are quite large in both GPD models. Although the value of £ is quite large, we see that the gluon 
contribution (i.e. the difference between the dashed and the solid curve) is by no means negligible. 

On Fig. [TUjwe show the DVCS observables relevant to the COMPASS experiment at CERN, namely (from left to 
right) the mixed charge-spin asymmetry, the mixed charge-spin difference and the mixed charge-spin sum defined in 
Eq. (59) of Ref. [H, at the kinematics £ = 0.05, Q 2 = 4 GeV 2 and -t = 0.2 GeV 2 . The upper part of Fig. [TO] uses 
the GK parametrization and the lower part the MSTW08 parametrizations with ler errors. We display here only the 
contribution from the GPD H . The lower value of the skewness £ = 0.05 allows to test a complementary regime with 
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Figure 4: The real part of the timelike Compton Form Factor H multiplied by 77, as a function of rj in the double distribution 
model based on Kroll-Goloskokov (upper left) and MSTW08 (upper right) parametrizations, for fi% — Q 2 — 4 GeV 2 and 
t — —0.1 GeV 2 . Below the ratios of the NLO correction to LO result of the corresponding models. 




0.2 n ■ — ■ — ■ — 1 - ■■! I 0.2 rr 




1.6 u ' ' J - 0.6 u 1 1 L - 

10- 3 10~ 2 10-' n 1 10~ 3 10- 2 10-' „ 1 



Figure 5: The imaginary part of the timelike Compton Form Factor W multiplied by rj, as a function of rj in the double 
distribution model based on Kroll-Goloskokov (upper left) and MSTW08 (upper right) parametrizations, for = Q 2 = 4 GeV 2 
and t — —0.1 GeV 2 . Below the ratios of the NLO correction to LO result of the corresponding models. 
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Figure 6: The ratio of the timelike to spacelike NLO corrections in the real (left) and imaginary (right) part of the Compton 
Form Factor H, as a function of £ in the double distribution model based on Kroll-Goloskokov (dashed) and MSTW08 (solid) 
parametrizations, for fip — Q 2 = 4 GeV 2 and t = —0.1 GeV 2 . For comparison timelike CFFs where calculated at Tj = £. 



0.0 




0.01 0.1 „ 1 



Figure 7: The real (left) and imaginary(right) parts of the TCS Compton Form Factor H multiplied by r/, as a function of rj 
in the double distribution model based on MSTW08 parametrization, for /i 2 ? = Q 2 = 4 GeV 2 and t = —0.1 GeV 2 . The dotted 
line shows the LO result and shaded bands around solid lines show the effect of a one sigma uncertainty of the input MSTW08 
fit to the full NLO result. 



respect to JLab measurements. Note the dramatic difference in the real part once gluon GPDs are in, inducing for 
instance a change of sign of the mixed charge-spin asymmetry and the mixed charge-spin difference for the MSTW 
case for instance. Note that the change is also huge in the case of GK parametrization but there is no sign change, 
indicating a significant model dependence. At any rate, NLO effects should be highly visible at COMPASS, which 
probes processes occuring at higher energy than JLab. 



2. Timelike Compton scattering 

Now we pass to predictions for the observables in the timelike counterpart of DVCS, namely TCS. On the left 
part of Fig. [TTJ we show the TCS contribution to the differential cross section as a function of the skewness r\ for 
Q 2 = M 2 = 4 GeV 2 , and t = -0.2 GeV 2 integrated over G (7r/4,37r/4) and over cf> G (0,2tt). We see that the 
inclusion of the NLO corrections is more important at small skewness. On the right panel of Fig. [TTJ we show 
that the Bethe-Heitler dominates the integrated cross-section for this kinematics. In consequence, more differential 
observables, as the azimuthal 4> dependence (with angles and </> defined in Ref. [7[ ) reveal in a better way the 
different contributions. Moreover simple <fi dependence of the interference term allows for an easy access to the real 
part of the CFFs, which as we observed on the Fig. [TJ are subject to the big NLO corrections. We indeed observe 
that effect on the Fig. [TU which shows the <f> dependence of the unpolarized differential cross sections for pure BH 
process, and with a LO and NLO corrections to the interference term. 

To quantify how big is the deviation from pure Bethe-Heitler process in the unpolarized cross section we calculate 
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Figure 8: The real (left) and imaginary (right) parts of the timelike Compton Form Factor T-L multiplied by 77 as a function of rj 
in the double distribution model based on KG parametrizations, for fj,% = Q 2 — 4 GeV 2 and t — —0.1 GeV 2 , at the Born order 
(dotted line), including the NLO quark corrections (dashed line) and including both quark and gluon NLO corrections (solid 
line) . The plots in the lower part show the corresponding ratios of the NLO correction to the LO results. 



(see Fig. [T3J) the ratio R, defined in Ref. [7| by 



2 ffdtp cos ip ds 



%) = ^7F-^' (31) 

J0 ®-<-P d Q,2 dtdip 

where S is a weighted cross section given by Eq. (43) of Ref. [7|. It is plotted on Fig. [13] as a function of the skewness ij 
for Q 2 = fj 2 = 4 GeV 2 , and t = —0.2 GeV 2 . In the leading twist the numerator is linear in the real part of the CFFs, 
and the denominator, for the kinematics we consider, is dominated by the Bethe - Heitler contribution. The inclusion 
of NLO corrections to the TCS amplitude is indeed dramatic for such an observable and includes also change of sign. 

Imaginary parts of the CFFs are accesible through observables making use of photon circular polarizations [7] . The 
photon beam circular polarization asymmetry 

A=4^, (32) 

is shown on the left part of Fig. [TH as a function of <fi for the kinematic variables relevant for JLab: Q 2 — 4 GeV 2 = 
lip, t = —0.1 GeV 2 and E 7 = 10 GeV (which corresponds to rj « 0.11). The same quantity is shown on the right 
panel of Fig. [L4l as a function of r\ for tfi — it/ 2 and Q 2 = 4 GeV 2 = /ip. The effect of the NLO corrections on that 
observable is rather large, ranging from 10% at the 77 = 0.1 (relevant for JLab) through 30% at 77 = 0.05 (relevant for 
COMPASS) up to 100% at very small rj's. 



6. CONCLUSION. 



Deeply Virtual Compton Scattering, both in its spacelike and timelike realizations, is the golden channel to extract 
GPDs from experimental observables. This extraction may be seen as a two step process: firstly, Compton form 
factors may be separated from a careful analysis of various differential cross sections and asymmetries. Secondly, 
convolutions of coefficient functions with model GPDs may be confronted to these CFFs. We have demonstrated 
here, in the case of medium energy kinematics which will be explored in the near future at JLab and COMPASS, 
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Figure 9: From left to right, the total DVCS cross section in pb/GeV 4 , the difference of cross sections for opposite lepton 
helicities in pb/GeV 4 , the corresponding asymmetry, all as a function of the usual <f> angle (in the Trento convention [33]) for 
E e = 11 GeV, fi 2 F = Q 2 = 4 GeV 2 and t = -0.2 GeV 2 . On the first line, the GPD H(x, £, t) is parametrized by the GK model, 
on the second line H(x,^,t) is parametrized in the double distribution model based on the MSTW08 parametrization. The 
contributions from other GPDs are not included. In all plots, the LO result is shown as the dotted line, the full NLO result by 
the solid line and the NLO result without the gluonic contribution as the dashed line. The Bethe-Heitler contribution appears 
as the dash-dotted line in the cross section plots (left part). 
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Figure 10: The DVCS observables for the COMPASS experiment, from left to right, mixed charge-spin asymmetry, mixed 
charge-spin difference and mixed charge-spin sum (in nb/GeV 4 ). The kinematical point is chosen as £ = 0.05, Q 2 = 4 GeV 2 , 
t — —0.2 GeV 2 . On the first line, the GPD H(x,£,t) is parametrized by the GK model, on the second line H(x,£,t) is 
parametrized in the double distribution model based on the MSTW08 parametrization. The contributions from other GPDs 
are not included. In all plots, the LO result is shown as the dotted line, the full NLO result by the solid line and the NLO 
result without the gluonic contribution as the dashed line. 
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Figure 11: left : LO (dotted) and NLO (solid) the TCS contribution to the cross-section as a function of r\ for Q 2 = /i 2 — 4 GeV 2 , 
and t = -0.2 GeV 2 integrated over 9 G (vr/4, 3vr/4) and over <j> G (0,2tt). Right : LO (dotted) and NLO (solid) TCS and 
Bethe- Heftier (dashed) contributions to the cross-section as a function of t for Q 2 = fi 2 GeV 2 integrated over 8 G (vr/4, 3n/4) 
and over G (0, 2tt) for £ 7 = 10 GeV (77 « 0.11). 
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Figure 12: The <f> dependence of the cross-section at E-y = 10 GeV, Q 2 = (j 2 = 4 GeV 2 , and t = —0.1 GeV 2 integrated over 
6 G (tt/4, 37t/4): pure Bethe-Heftler contribution (dashed), Bethe-Heitler plus interference contribution at LO (dotted) and 
NLO (solid). 



that the inclusion of NLO corrections to the coefficient function was an important issue, and that the difference of 
these corrections between the spacelike and timelike regimes was so sizeable that they can be promoted to the status 
of direct tests of the QCD understanding of the reactions. 

Let us stress again a feature that was largely overlooked in previous studies, namely the importance of gluon 
contributions to the DVCS amplitude, even when the skewness variable £ is in the so-called valence region. This is 
not a real surprise when one recalls that gluons (in terms of distribution functions) are by no means restricted to the 
very low x region and that gluon CFFs at a given £ value also depend on gluon PDFs at lower values of x. This 
effect is particularly big when one considers the real part of CFFs in the timelike case. This promotes the observables 
related to this quantity as sensitive probes of the 3-dimensional gluon content of the nucleon. 

We did not extend our study to the very high energy regime which will be explored (both for spacelike and timelike 
Compton scattering) at future electron-ion colliders, nor to the case of ultraperipheral collisions at present hadron 
colliders. This will be addressed in separate studies. We did not discuss the rich scope of factorization scale dependence 
issues, which deserves special attention. Contrarily to the case of inclusive reactions where various strategies have 
been built to optimize the factorization scale, it has been shown in Ref. (34j that it was quite impossible to find a 
recipe to minimize higher order corrections both for the real and the imaginary part of a CFF. Moreover we find it 
difficult to advocate a different choice of scale for the two cases of timelike and spacelike reactions. 
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Figure 13: Ratio R defined by Eq. (J3TJ as a function of 77, for Q 2 = fi 2 F = 4 GeV 2 and t = —0.1 GeV 2 . The dotted line 
represents LO contribution and the solid line represents NLO result. 




Figure 14: (Left) Photon beam circular polarization asymmetry as a function of cj>, for t = —0.1 GeV 2 , Q 2 = /i 2 = 4 GeV 2 , 
integrated over 9 £ (tt/4, 37r/4) and for _E 7 = 10 GeV (77 ~ 0.11). (Right) The rj dependence of the photon beam circular 
polarization asymmetry for Q 2 — fi 2 = 4 GeV 2 , and t = —0.2 GeV 2 integrated over 9 € (tt/4, 37r/4). The LO result is shown 
as the dotted line, the full NLO result by the solid line. 
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